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Abstract 

A detailed consideration is given to the translation-invariant theory of Tulub polaron 
constructed without the use of Pekar ansatz. A fundamental result of the theory is that 
the value of the polaron energy is lower than that obtained on the basis of Pekar ansatz 
which was considered as an asymptotically exact solution in the strong coupling limit. 
In the case of bipolarons the theory yields the best values of the coupling energy and 
critical parameters of their stability. Numerous physical consequences of the existence of 
translation-invariant polarons and bipolarons are discussed. 
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1 Introduction. Pekar ansatz. 


As is known, polaron theory was among the first to describe the interaction between a particle 
and a quantum held. Various aspects of the polaron theory are presented in numerous reviews 
and books ra-nm- Being non-relativistic, the theory does not contain any divergencies and 
for more than sixty years has been a testing area for approbation of various methods of the 
quantum held theory. Though no exact solution of the polaron problem has been found up to 
now, it has been believed that the properties of the ground state are known in detail. This 
primarily refers to the limit cases of weak and strong coupling. A solution in the weak coupling 
limit was given by Froehlich HU, and that in the strong coupling one was found by Pekar 
cun. By now rather an exact solution has been obtained for the energy of the polaron 
ground state in the weak coupling limit [13J, na: 

E 0 = -(a + 0.0159196220a 2 + 0.000806070048a 3 -) , (1.1) 

where hu 0 is the energy of an optical phonon, a is the constant of electron-phonon coupling. 

A solution of the problem in the opposite strong coupling limit was given by Pekar on the 
assumption that the wave function T of the electron + held system has the form: 

T(r, q u ..., q u ...) = ^(r)$(g i,..., q % , ...) , (1.2) 

where r ) is the electrons wave function depending only on the electron coordinates, $ - 
is the wave function of the held depending only on the held coordinates. Pekar himself [Ij 
considered ansatz (H3 ) to be an approximate solution. In the pioneer works by Bogolyubov 
and Tyablikov |15].[T6] it was shown that in a consistent translation-invariant theory the use 
of ansatz (II.2p (for decomposed coordinates introduced in [15]. [15] ) gives the same results for 
the polaron ground state energy as the semiclassical Pekar theory does am With the use 
of (II . 2[) the value of the ground state energy has been found to a high precision. According to 
HE, IE] it is equal to:: 

E = (-0.108513a 2 - 2.836) huj 0 ■ (1.3) 

The idea that Pekar ansatz (11.21) is an exact solution of the strong coupling polaron problem 
was completely established after publication of [P9] where asymptotics (jl ,3j) was strictly proved 
by path integral method, i.e. without the use of ansatz (j 1.2 j) (see also review [2U]). 
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Before the publication of paper [TU] many attempts were made to improve the strong cou¬ 
pling theory |2T]- |E6] • The reason why Pekar ansatz caused the feeling of disappointment was 
translation invariance of the initial polaron Hamiltonian. When ansatz for the wave function 
"0(r) (II.2p is used the wave equation has a localized solution. The electron is localized in a 
potential polarization well induced by it. In other words, the solution obtained does not possess 
the symmetry of the initial Hamiltonian. Self-trapping of the electron in the localized potential 
well leads to a spontaneous breaking of the systems symmetry. Attempts to restore the initial 
symmetry were based on the use of degeneration of the system with broken symmetry. Since 
in a homogeneous and isotropic medium nothing should depend on the position of the polaron 
well center r 0 , one can ”spread” the initially localized solution over all the positions of the 
polaron potential well by choosing the wave function in the form of a linear combination in all 
the positions of the well. 

In the most consistent form this program was carried out in (24] . With this end in view 
for the wave function which is an eigen function of the total momentum, the authors used a 
superposition of plain waves corresponding to the total momentum multiplied by wave functions 
obtained from (11.21) to which a translations operator is applied. In other words, they took an 
appropriate superposition with respect to all the positions of the polaron well r 0 . The main 
result of paper m is that calculation of the polaron ground state energy with such a delocalized 
function yields the same value as calculations with localized function (1 1.2 [) do. The authors 
of [23J also reproduced the value of the polaron mass which was earlier obtained by Landau 
and Pekar ra on the assumption that polaron moves in medium in the localized state (11.21) . 
The results derived in m were an important step in resolving the contradiction between the 
requirement that the translation-invariant wave function be delocalized while the wave function 
of the self-trapped state be localized. 

Notwithstanding the success achieved with this approach it cannot be considered fully ade¬ 
quate since it has quite a few inconsistencies. They follow from the very nature of the semiclas- 
sical description used. Indeed, the superposition constructed in m on the one hand determines 
the polaron delocalized state, but on the other hand, without changing this state, one can mea¬ 
sure its position and find out a localized polaron well with an electron localized in it. The 
reason of this paradox is a classical character of the polaron well in the strong coupling limit 
and, as a consequence, commutation of the total momentum operator with the position of the 
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polaron well r 0 . □ 

To remedy this defect some approaches were suggested in which the quantity r 0 (r, qi ,..., q i: ...) 
which is not actually an additional degree of freedom was considered to be that though with 
some additional constraints. Discussion of these challenges associated with solution of the prob¬ 
lem of introducing collective coordinates is given in review [301]. Since the results obtained by 
introducing collective coordinates into the polaron theory are polemical it seems appropriate 
to describe strict results of the translation-invariant theory without recourse to the concept 
of collective coordinates. The aim of this review is to present an approach used in the strong 
coupling limit which does not use Pekar ansatz. 

A solution possessing these properties in the case of a strong coupling polaron was originally 
found by Tulub [3TJ, [32] • During nearly half a century the result obtained in [3T],[32j was n °f 
recognized by specialists working in the field of polaron theory. The reason why the importance 
of that result was not appreciated was an improper choice of the probe wave function in [32j 
to estimate the ground state. As a result the ground state energy was found in [32] to be 
Eq = —0.105a 1 2 /kao which is larger than in (11.31) . An appropriate choice of the wave function 
has been made quite recently in paper [33]. This has yielded a lower than in (1 1.3 [) value of 
the polaron ground state energy equal to E 0 = —0.125720 a 2 hu 0 . Hence, actually we have 
to do with inapplicability of adiabatic approximation in the case of a polaron, though it is 
fundamental for solid state physics. 

In this review we present the main points of the translation-invariant polaron (Tl-polaron) 
theory. In §2 Heisenberg canonical transformation is used to introduce a coordinate-free Pekar 
Froehlich Hamiltonian which forms the basis for translation-invariant description. With the 
use of Lee-Low-Pines wave function the weak coupling limit is reproduced. 

In §3 we present the general translation-invariant Tulub theory valid for any values of 
the electron-phonon coupling constant. In view of primary importance of Tulub approach the 
material is presented in greater detail than in the original paper so that a reader may reproduce 
and check the results obtained if desired. A general expression is obtained for the Tl-polaron 

1 At the rise of quantum mechanics, the founders of the science were fully aware of the difficulties arising here. 

Thus, for example, in m Bethe notices that for a proper quantum-mechanical description of an interaction 
between a field and particles, quantizing of the field is required, i.e. quantum theory of the field: ” The 
fact is that, when quantizing mechanical parameters (coordinates and momenta) one should also quantize the 
associated fields. Otherwise, as Bohr and Rosenfeld showed j2^, an imaginary experiment can be suggested 
which consists in simultaneous measurement of the coordinate and the momentum of a particle from examination 
of the field induced by it. This contravenes Heisenberg’s Uncertainty Principle.” 
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energy for an arbitrary coupling strength. 

In §4 a limit case of weak coupling is considered. It is shown that the general expression 
for the polaron energy derived in the previous section reproduces the weak coupling limit to a 
high precision. 

In §5 a limit case of strong coupling is dealt with. It is shown that in this limit case the 
polaron ground state energy has a lower value as compared to that obtained on the basis of 
Pekar ansatz. 

In §6 the translation-invariant polaron theory is generalized to the case of a bipolaron. 
The bipolaron ground state energy obtained here is much lower than that yielded by the best 
variational calculations with the use of Pekar ansatz. The calculation results are used to explain 
the high-temperature superconductivity by using Tl-bipolarons. 

In §7 the energy of Tl-polaron and bipolaron is derived in an alternative way. The approach 
presented there enables one to get an explicit form of wave functions for a polaron and bipolaron. 
The results obtained suggest that for all the values of the electron-phonon coupling constant 
the wave functions derived describe delocalized states. Hence it is shown that in the strong 
coupling limit Pekar ansatz is not fulfilled and no transition to the self-trapped (i.e. localized) 
state with broken symmetry takes place. 

The results obtained radically change the visions of polarons and bipolarons and in the 
general case cast some doubt on the concept of self-localized states in condensed systems. 
Recently an open discussion on the completeness of Tulub theory has taken place. The main 
points and results of the discussion are presented in §8. 

In §9 we discuss some problems and consequences emerging from the existence of Tl-polarons 
and bipolarons. Of special practical interest are results concerned with their superconducting 
properties. 

§10 is devoted to discussion of some fundamental problems of the strong coupling theory 
which are still to be solved. 

In Appendices 1-3 we give proofs of some important statements basic to the approach under 
consideration. 
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2 Coordinate-free Hamiltonian. Weak coupling. 


Let us proceed from Pekar-Froehlich Hamiltonian: 

h 2 


H = A r + J2 '4 ( a k e‘ kr + 


( 2 . 1 ) 


where a^, a k are operators of the birth and annihilation of the field quanta with energy hui® = 
huio, m is the electron effective mass, 14 is a function of the wave vector k. 

Interest in the study of this Hamiltonian is also provoked by the fact that, as distinct 
from many model Hamiltonians considered in the condensed matter theory, Pekar-Froehlich 
Hamiltonian (12. lj) in the limit of long waves asymptotically exactly describes the behavior of a 
non-relativistic electron in a continuous polar medium. 

Electron coordinates can be excluded from (12.ip via Heisenberg transformation [34] : 

Si = exp (p - y, hka^a^j r j , (2.2) 

where P is the total momentum of the system. Application of S\ to the field operators yields: 

S^a k S x = a k e~ ikr , Af'a+A, = a+e ikr . 

Accordingly the transformed Hamiltonian H = S^HSi takes on the form: 

H = ^-(p-^hka^a^] +5^14 (a k + a^) + ^2nw° k a^a k . (2.3) 

m \ k J k k 

Since Hamiltonian (12. 3 p does not contain electron coordinates, it is obvious that solution of the 
polaron problem obtained on the basis of (12.3p is translation-invariant. Lee, Low and Pines 
[ 35] studied the ground state (12.31) with the probe wave function \^>)llp '■ 


\^)llp — S 2 |0) j 


(2.4) 


where: 


S 2 = 


exp |- at)! , 


(2.5) 


f k are variational parameters having the meaning of the value of displacement of the held 
oscillators from their equilibrium positions, |0) is the vacuum wave function. The quantity f k 
in S *2 (12.5p is determined by minimization of energy E = (01/S^ -1 10), which for P = 0 yields: 


E = 2j2fkV k + 


2m 


w, 


£^a 2 + £^a 2 . 


2 771 

k k 


( 2 . 6 ) 
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(2.7) 


ft 


In the case of an ionic crystal: 


Vk 

hu ° + K 2 k 2 /2m ' 


e l2TThuJo 

Vk = - k y— 


hujQ /47ra\^ 2 f2mu>o\ 1 ^ 2 

ku V2 \ V~) ’ u = V h ) 


1 e 2 u 

2 hjj 0 e ’ 


p- 1 = p- 1 -f- 1 
t oo fc 0 


( 2 . 8 ) 


where e is an electron charge, and £ 0 are high-frequency and static dielectric permittivities, 
a is a constant of electron-phonon coupling. With substitution of (12.8)) into (12.61) . (12.71) the 
ground state energy becomes E = —ahuj o, which is the energy of a weak coupling polaron in 
the first order with respect to a. 

A solution of the problem of transition to the strong coupling case in coordinate-free Hamil¬ 
tonian (12. 3 p was found on the basis of the general translation-invariant theory constructed in 
Tulubs work [52] , The main points of this theory are given in the next section. 


3 Coordinate-free Hamiltonian. General case. 


To construct the general translation-invariant theory in works of [31] , [32] was used a canonical 
transformation of Hamiltonian ()2.3[) with the use of operator S 2 ()2.5|) which leads to a shift of 
the field operators: 

S 2 1 a k S 2 = a k +f k , S 2 1 alS 2 = a++ f k . (3.1) 

The resultant Hamiltonian H = S,J 1 H S 2 has the form: 

ti = H 0 + H u (3.2) 


where: 


Hn = 


P 2 
2 m 


2 £ua+£ (ms- 


hkP 


m 


fl + 


2m 


£*71 +«o, 


(3.3) 


Ho = £ hwta^at + ^— £ kk'ftfk' {dkdf + a)p£ + a*(if + a^at) , (3.4) 


hjjju = fvjJi 


kM' 


kkP h 2 k 2 hk 


m 


+ 


+-£»-/y 

m ^^ 


2m m 


(3.5) 


k' 


Hamiltonian Hi contains terms linear, triple and quadruple in the birth and annihilation op¬ 
erators. With an appropriate choice of the wave function diagonalizing quadratic form (13. 4 p 
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mathematical expectation Hi becomes zero (Appendix 1). In what follows we believe that 
h—1, ojq = 1, m = 1. To transform Hq to a diagonal form we put: 




Qk = -yL= ( a k + flfc) , Pk = “ a k) > z k = kfkV2uT k - 

With the use of (13. 6 p expression (j3.4[) is written as: 

n v = \j2( p kP k+uj fak ( ik)+ ^{j2* k q k ) 

k \ k ) k 

This yields (13. 7 j) the following motion equation for operator q k : 

q k + co 2 k q k = -z k ^ z k /q k i . 

k' 

Let us search for a solution of system (13.8p in the form: 

Qk(t) = 22 ^kk'Ck'(t) , Ck(t) = . 

k' 

As a result we express matrix Q kk j as: 

(vu - ul)Stw = Z t Zy,a 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


‘‘k"k' 


(3.10) 


Let us consider determinant of this system which is derived by replacing the eigenvalues v\ in 
(I3.10p with the quantity s which can differ from u 2 . The determinant of this system will be 


det | (s - ul)5 kkl - z k z k ' \ = - i%). 

k 

On the other hand, according to [36: 


(3.11) 


za 


z 


k' 


det | (s - ul)5 kk > - z k z k /1 = TT(s - u\) 1 - - S' 

k \ 6 k' s 

It is convenient to introduce the quantity A(s) : 

a (s)=- u i)/ nc- - ^ • 

k k 

With the use of (13.lip and (I3.12p A(s) is expressed as: 

A w= fi-sE 


(3.12) 


(3.13) 


3 “ s - ui‘ k 


(3.14) 


2 In Wentzel’s work z k is not a vector function, but a scalar one, therefore a ’’cube” in (13.121) is lacking. 
Generalization to the vector case is given in m- 


















From (13.lip . (I3.12p follows that the frequencies renormalized by interaction are determined 
by a solution to the equation: 

AK 2 ) = 0 . (3.15) 

The change in the systems energy A E caused by the electron-held interaction is equal to: 

A E = ^'^2 / (v k -u k ) ■ (3.16) 

k 

To express the quantity A E via A(s) we use Wentzel approach |36]. Following [36] we write 
down the identity equation: 






1 

2ni 


dsf(s)— lnA(s) 


1 

2 m 


dsf\s) In A(s), (3.17) 


c c 

where integration is carried out over the contour presented in Fig.l 



Fig. 1. Contour C. 


Taking f(s) = y/s we get: 


A E = 


ds 


In A(s) 


c 


Turning in (I3.14p from summing up to integration with the use of the relation: 

1 


E 


(2 7T) 


d k 


in a continuous case, using for expression (I3.14p for A(s) we obtain: 

2 f k 2 flu 2 


A (s) = D 6 (s), D(s) = 1 — 


I k^k 73 


3(27r) 3 ,/ s — (jj\ 


d 6 k. 


(3.18) 


(3.19) 
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As a result the total energy of the electron is: 


E = AE + 2^ V k f k + ^ f 2 k u° k . (3.20) 

k k 

The results obtained here are general and valid for various polaron models (i.e. any functions V k 
and oo k ). Below we consider limit cases of weak and strong coupling which follow from general 
expression (13.201) on the assumption that P = 0. 

Notice that for P/ 0, according to [31] , expression (13.20j) takes the form: 


P 


e =^ +ae( ~ p '> +2 E, v ^ + E\^ 


hkP 


771 




1 1 . 


A E(P) = 


3 

—— / -^hiTTzr(s), 

8m j 1 ; ’ 


D\S) = 1 - £ 


2 ’ 

* s ~ w “ 


where z\ - i -th component of the vector z k . Functions f k , ui k and z k are depending on \k\ and 
on ( kP ). 


4 Weak coupling limit in Tulub theory. 

Quantities f k in the expression for the total energy E should be found from the minimum 
condition: 5E/5f k = 0 which yields the following integral equation for f k : 


fk = -V k /(1 + k-/2/j, k ), n k = 


_ i Uk f ds 


2m J y/s (s — ojI)D(s) ’ 


(4.1) 


c 


In the case of weak coupling a —> 0 and equations (14. Tj) can be solved with the use of pertur 
bation theory. In a first approximation as a —> 0, D{s) = 1 and /i^ 1 is equal to: 

c o k f ds 1 


k-k ~ 


2iri J y/s (s - c o k ) 


= 1. 


(4.2) 


c 


Accordingly, f k from (14.11) is written as: 


fk = -V k /( 1 + k 2 /2 ). 


(4.3) 
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The quantity A E involved in the total energy takes on the form: 

AE — —— J> in D{s) , l„Z>( s ) = --i- f^±cPk. (4.4) 

8m J ^ y y J 3(2vr) 3 J s-co 2 k 

c 

With the use of (14. 3 p integrals involved in (j4.4j) are found to be: AE = ( a/2)hu 0 . Having 
calculated the rest of the terms involved in expression (13.20[) we get the first term of the 
expansion of polaron total energy in the coupling constant a : E — —ahujQ. 

In papers |3T], [37], [38] a general scheme of calculating the higher terms of expansion in a 
was developed. In particular, the eigen energy and effective mass were found to be [38] : 

E = — (a + 0.01592a 2 )/kno , 


a 


m* = (1 + - + 0.02362a ) m . 


(4.5) 


Hence within the accuracy of the terms 0(a 3 ) the polaron energy expression calculated within 
Tulub approach with the use of perturbation theory coincides with exact result (11.11) (see §7). 


5 Strong coupling. 


The case of strong coupling is much more complicated. To reveal the character of the solution in 
the strong coupling region let us start with considering the analytical properties of the function 
D (s) in the form: 

OO 

s — 1 f k A flukdk 


D(s) = D(l) + 


3^ J K-1)K- S )’ 


(5.1) 


where D( 1) is the value of D(s) for s — 1: 


77(1) — 1 + Q = 1 + 


OO 

1 f k4 fk^k 


3n 2 J col — 1 
o 


dk . 


(5.2) 


From (I3.19|) also follows that: 


£>(») = !- 




3n 2 J s — col 


dk. 


(5.3) 


Function D (s), being a function of a complex variable s, has the following properties: 1) 
D(s) has a crosscut along the real axis from s — 1 to oo and has no other peculiarities; 2) 
D*(s ) = D(s*)] 3) as s —» oo sD(s) increases not slower than s. These properties enable us to 
present the function [(s — l)Z7(s)] _1 in the form (Appendix 2): 

1 1 ' ds ' (5.4) 


(s — l).D(s) 2ni J (s' — s)(s' — l)D(s') ’ 

C~\~P 
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where contour C + p is shown in Fig. 2: 



Fig. 2. Contour C + p. 


The integrand function in (j5.4[i has a pole at s' = 1 and a section from s' — 1 to s' — oo. 
Having performed integration in (j5.4jl along the upper and bottom sides of the crosscut we get 
the following integral equation for Zl _1 (s) : 


1 s- 1 

+ 


k^flujkdk 


2M2 • 


D(s) 1 + Q 3^ J {s-ul){ul-l)\D(ul)\ 

o 

With the use of integration by parts expression (13.18j) for A E can be written as: 


A E = 


dkk A flu k 


2vr 2 


vs 


27 n J (s — tu|) 2 D(s) 
c 


ds 


From (15.51) . (15.61) we have: 


A E = 


k 4 P k dk + 1 


2t r 2 J 2(1 + Q) 12t r 4 



0 0 


k A pP 4 fpU P {u k u P + c o k (u: k + uj p ) + 1 ) 
(u k +up) 2 (u 2 - l)\D(ul )\ 2 


Equation ji k (14.1)) . according to (15. 5 F can be presented in the form: 


-i 1 1 

d k ~ ^ + 


P A fp(u k ujp + 1 )dp 


1 + Q 3n 2 J (uj 2 -l)(uj k +ujp)\D(uj 2 )\ 2 


(5.5) 


(5.6) 


dpdk. (5.7) 


(5.8) 


Equations (14.ip . (15.8[) for finding f k as well as expressions (I3.20[) . (15.71) for calculating polaron 
energy are very complicated and their exact solution can hardly be obtained. To calculate 
approximately the energy E given by (13.201) , (15.71) in [22] a direct variational principle was 
used. For the probe function, the author used Gaussian function of the form: 


fk = ~V k exp(—k 2 /2a 2 ), 


(5.9) 
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where a is a variable parameter, besides, as can be seen in the case of strong coupling, a 1. 
Substitution of (15. 9 p into (j3. 19j) yields for real and imaginary parts of D(s) (see Appendix 3): 

Re D(ul) = 1 + A v(y ), ImD(^) = k 3 fl/ 6tt , 


v(y) = 1 — ye y / e l dt — ye y / e f dt , 


(5.10) 


A = 4cra/3\/27r, y = k/a. 

In the limit of strong coupling (ex 1) the expression for energy E given by (13. 201) with the 
use of (I5.7j) takes on the form: 


s -ii- 




A ,= 


1 + (? < A 


f e-y\l-n(y))dy 
tt J (1/A + v(y)) 2 + Try 2 e~ 2y2 /A ’ 


(5.11) 

(5.12) 


Sl(y) = 2y‘{ (1 + 2y V”' J e^dt-y 

y 

As A —> exo, integral (I5.12p has maximum for y 4 = 3A/4, if the function / fc is chosen in the form 
(15.91) . however if the actual boundedness of the region of integration with respect to y is taken 
into account, this peculiarity does not take place (see §8). 

When calculating (j5.12j) in paper [32] Tulub assumed that in the strong coupling limit 
1/A = 0. As a result of numerical integration g(0) was found to be q(0) = 5.75 whence, varying 
energy E (15.lip with respect to a we get: 


E = —0.105cx 2 ho;i 


o • 


(5.13) 


Comparison of (15.13p with (11.31) shows that the value of E obtained for a —» oo lies higher than 
the exact value in Pekar theory (11.31) . For this reason, until quite recently it was believed that 
Tulub theory as applied to a polaron does not give any new results. 

The situation changed radically after publication of [33]. There it was shown that the choice 
of the wave function for minimizing energy (I3.20p in the form (15.9p is not optimal since it does 
not satisfy viria.l relations. As is shown in [33], an appropriate function / fc should contain the 
multiplier \/2 outside the exponent in expression (15.9p . Minimization of energy (13.201) with the 
optimal probe function yields: 

E = —0.125720cx 2 hcn 0 . (5.14) 
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Result (15. 14j) is fundamental. Above all it means that Pekar ansatz does not give an exact 
solution. Though result (I5.14p refers to a particular case of PekarFroehlich Hamiltonian with 
Vk given by (12.8)) . the conceptual conclusion should be valid for all types of self-localized states. 
Of special interest is to consider the case of bipolarons since they can play an important role 
in superconductivity. 

6 Tl-bipolarons. 

Great attention given to polaron problem in recent times is associated with attempts to explain 
the superconductivity phenomenon relying on the mechanism of Bose condensation of bipolaron 
gas. In this connection the study of conditions under which the bipolaron states are stable is 
of paramount importance. The theory of large-radius bipolarons which are now considered to 
be the best candidates for the role of charged bosons forming Bose-Einstein condensate with 
pairing in real space is considered in detail in reviews 0. 0. H- 

The study of the process of formation of a stable two-electron state in a crystal, or a 
bipolaron, generally implies finding a pairwise interaction between two polarons as a function 
of a distance between them [8]. For a large-radius bipolaron, the region of its existence is 
bordered on the part of the coupling constant a by rather a large value of a c below which the 
polaron bound state does not exist. In view of the requirement that a c be large, which may not 
be met in high-temperature superconductors, some researchers investigated the contribution of 
other types of interactions and coupling symmetries [3D], [31] • 

In what follows we will deal with only electron-phonon PekarFroehlich interaction since the 
approach under consideration can be generalized to other types of interactions as well. It seems 
all the more actual in view of the fact that in recent times some reasoned arguments have been 
obtained testifying that electron-phonon interaction in high-temperature superconductors is 
strong [32]- [S3] - There are also arguments in favor of the fact that owing to weak screening of 
high-frequency optical phonons, the electron-phonon interaction in high-temperature supercon¬ 
ductors is more adequately described not in the framework of a contact interaction of Holstein 
polaron model 0 a, but in terms of a long-range interaction of Froehlich type [46]. 

Before the publication of W\-W\ the lowest values of the energy of bipolaron states deter¬ 
mined by electron-phonon interaction were obtained in [I5QH5Z] for a < 8 and in [52j-[55] for 
a > 8. Attempts to find a translation-invariant solution of the bipolaron problem by varia- 
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tional methods using direct variation of the wave function of the two-electron system [32], [5jB], 
|57] yielded higher values of the bipolaron ground state energy as compared to those obtained 
with the use of a wave function lacking translation invariance [51], [52], [55], [58] - In this sec¬ 
tion we present the results obtained in [T7]-[42]J for a bipolaron within the translation-invariant 
approach. 

Let us proceed from Pekar-Frochlich Hamiltonian for a bipolaron [8] : 


H = -^A ri - ^A r2 + hiu 0 k a+a k + U (|f=i - r 2 |) + ^ (V k e ikr 'a k + V k e ikr *a k + H.C .) , 

k k 

(6- 1 ) 

U(\r l -r 2 \)= f 

£oo|n - r 2 1 

where f\ and r 2 are coordinates of the first and second electrons, respectively, quantity U 
describes Coulomb repulsion between the electrons. 

In the mass center system Hamiltonian (16.Tj) takes on the form: 


f ) 2 h 2 kr / -- \ 

H = ~^A R - — A r + U(\r\) + ^ fuj° k a+a k + cos — (a k e ikR + H.C.) , (6.2) 

" e k k 


2 M P 


R = (?i + r* 2 )/2 , r — f\ — r 2 , M e = 2m , /i e = m/2 . 

In what follows we will believe that h = 1, u k = 1, M e = 1 (accordingly /x e =1/4). 

The coordinates of the mass center R can be excluded from Hamiltonian (16.2[) via Heisenberg 
canonical transformation: 



__ _ _ Lj? ^ \ 

H = S^HSi = -2A r + C/(|r|) + ^a^a fe + ^2I4cosy(a i; + a^) + - \ S^ka//a k ) . (6.3) 

k k Vi k J 

From formula (16.3ft follows that the exact solution of the bipolaron problem is determined by the 
wave function -0(r) containing only relative coordinates r and, therefore, possessing translation 
invariance. 

Averaging of H over r i/{r) leads to the Hamiltonian H\ 

H — — I ka k a k J + a k a k + V k (a k + a, k ) + T + U , (6-4) 

V k ) k k 


V k = 2V k {i)\cos—\4,), U — (ip\U (r)\4>), T=- 2<^|A r |V>) ■ 
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Hamiltonian (16.4j) differs from Hamiltonian (12,3j) in that the quantity 14 in (12. 3 p is replaced by 
V k and constants T and U are added. Therefore, repeating the derivation performed in §3 we 
express the bipolaron energy E bp as: 


E bp = AE + 2j2Vkfk + J2ti + f + tJ i ( 6 - 5 ) 

k k 

where A E is given by expression (I5.7[i . From (|6. 5[i we can get expressions for the bipolaron 
energy varying E bp with respect to f k and ip. Since the equations obtained in this way are 
difficult to solve, for actual determining of the bipolaron energy we use a direct variational 
approach, assuming [39] : 


fk = —NVk exp(— /c 2 /2/t) , 
ip(r) = (2/7r£ 2 ) 3 / 4 exp(-r 2 /f), 


( 6 . 6 ) 


where N, y, i are variational parameters. For N = 1, expression ( 16 . 61 ) reproduces the results 
of work m, and for N — 1 and y —* oo , those of work [381J . 

Having substituted ( 16 . 6 |) into the expression for the total energy and then minimized the 
expression obtained with respect to parameter N we write E as: 


E(x, y-rj) = $(x, y, y)a 2 , (6.7) 

6 20.25 16 Jx 2 + 16 y 4^2/^ 

Hx, y ; v) = — + 2 - r( 2 7 + —r -r • 

x A x- + lby \/vr(a; z + 8 y) x(l — y) 

Here x, y are variational parameters: x = £a, y = a 2 //i, y = Soo/e^Eet us write & m i n for the 
minimum of function $ of parameters x and y. Fig.3 shows the dependence of &mi n on the 
parameter y. Fig.4 demonstrates the dependence of x min , y min on the parameter y. 



Fig. 3. Graph $ m in(y)- 
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Fig. 4. Graphs x min (r]), y min (rj),. 


Fig.3 suggests that E min (rj = 0) = —0.440636a 2 yields the lowest estimation of the bipolaron 
ground state energy as compared to all those obtained earlier by variational method. Horizontal 
lines in Fig.3 correspond to the energies: E\ = —0.217a 2 and E 2 = —0.2515a 2 , where E\ = 
2 E pi , E pi is the Pekar polaron ground state energy (11.31) : E 2 = 2 E p2 where E p2 is the ground 
state energy of a translation-invariant polaron (15.141) . Intersection of these lines with the curve 
Emin(jl) yields the critical values of the parameters y — r} ci — 0.3325 and rj = r] C2 = 0.289. 
For 77 > r ) C2 , the bipolaron decays into two translation-invariant polarons, for 77 > rj ci , it breaks 
down into Pekar polarons. The values of minimizing parameters x m i n and y m in for these values 
of 77 are x min (ti) = 5.87561, 7/ min (0) = 2.58537, a; min (0.289) = 8.16266, 7/ min (0.289) = 3.68098, 
x min (0.3325) = 8.88739, y mm (0.3325) = 4.03682. 

The critical value of the electron-phonon coupling constant a, determined from comparison 
of the energy expressions in the weak coupling limit (doubled energy of a weak coupling polaron: 
E = —2ahu 0 ) and in the strong coupling limit (E = —0.440636a 2 fo<;o)> at which the translation- 
invariant bipolaron is formed is equal to a c ~ 4.54, being the lowest estimate obtained by 
variational method. It should be emphasized that this value is conventional. Hamiltonian (16. 4 p 
coincides in structure with one-electron Hamiltonian (12.31) . therefore, as in the case of a polaron, 
the bipolaron energy, by ra. is an analytical function of a. For this reason at the point a = a c 
the bipolaron energy does not have any peculiarities and the bipolaron state exists over the 
whole range of a and 77 variation: 0 < a < 00 ; 0 < 77 < 1 — 1/2 y/2, for which E < 0. To solve 
the problem of the existence of a c value at which the bipolaron state can decay into individual 
polarons, one should perform calculations for the case of intermediate coupling. In particular, a 
scenario is possible when the bipolaron energy for some values of 77 will be lower than the energy 
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of two individual polarons for all values of a, i.e. the bipolaron state exists always. Notice, 
that for the derived state of the translation-invariant bipolaron, the virial theorem holds to a 
high precision. 

The problems of arising high-temperature superconductivity (HTSC) and explaning this 
phenomenon by formation of bipolaron states were dealt with in a number of papers and reviews 
0 . [59], [60] . In these works the existence of HTSC is explained by Bose condensation of 

a bipolaron gas. The temperature of Bose condensation T 0 = 3.3lh 2 n^ 3 /k B rrib p , which is 
believed to be equal to the critical temperature of a superconducting transition T c for m bp ~ 
10 m, depending on the bipolarons concentration no, varies in a wide range from Tq = 3 K at 
n = 10 18 cm -3 to T 0 pz 300/L at n ~ 10 21 cm -3 . In the latter case the bipolarons concentration 
is so high that for a bipolaron gas as well as for Cooper pairs, the compound character of a 
bipolaron when it ceases to behave like an individual particle should show up. In the case 
of still higher concentrations a bipolaron should decay into two polarons. According to (16.61) 
the characteristic size of a bipolaron state is equal to l and in dimension units is written as: 
i C orr = K 2 £x{jj)/me 2 . Here d corr has the meaning of a correlation length. The dependence x(r/) 
is given by Fig.4. Fig.4 suggests that over the whole range of 77 variation where the bipolaron 
state is stable the value of x changes only slightly: from x{rj = 0) ~ 6 to x(r] = 0.289) ~ 8 . 
Hence, even for 77 = 7 / c , the critical value of the concentration at which the bipolarons multipiece 
character is noticeable is of the order of n c = 10 21 cm“ 3 . This result testifies that a bipolaron 
mechanism of HTSC can occur in copper oxides. 


7 Ground state functional. Tulub ansatz. 


To diagonalize quadratic form (13. 4 p we can use Bogolyubov-Tyablikov transformation [ 6 Tj . Let 
us write a k for the operators of physical particles in which TL 0 is a diagonal operator. 

Let us diagonalize the quadratic form with the use of the transformation: 


a k = M lkk ,a k > + M 2 kk' a t' > 

k' k' 

a k =Y1 M *kk' a k' + M 2 kk ' a k' , 

k' k' 

so that the equalities: 


(7.1) 


[ofc, Ofe/] — [cXfcj Ctfc'] — fikk’ , [^ 0 ) a k ] “ UkOlk ■ 


(7.2) 
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be fulfilled. 


(7.3) 


Relying on the properties of unitary transformation (17.11) we have: 

m 2 m+ = m*m 2 t , 

(M+)- 1 = M x - m;(m*)- 1 m 2 . 

With the use of (17.31) transformation of the operators reciprocal to (17.11) takes on the form: 

a k — ^ M lkk ' a k' ~ ^2 M 2 kk ,a t' ) 

(7.4) 


k' 


k' 


at = 


Mikk'^t' ~ ^2 M2kk'dk' ■ 


According to 


k' k' 

matrices Mi and M 2 are: 


(M h2 ) kk ' = ~(u k u k >) 1/2 (^fc ± w*/) 


S(k - k') + (kk’)f k f k . 


c 2‘{u k Uk i ) 1 / 2 


(ul,-ul±ie)D ± (ul) 


(7.5) 


^±K) — i + 


& 4 u k 


37r 2 y - aE =p ie 


— dk, 


where the superscript sign in the right-hand side of (17.51) refers to Mi and the subscript sign 
to M 2 . As a result of diagonalization quadratic form (13.4[) changes to: 


Hq = A E + ^2 u kat 


a k . 


(7.6) 


Functional of the ground state A 0 |0) is chosen from the condition: 


a fc A 0 |0) = 0 


(7.7) 


The explicit form of functional Ao is conveniently derived if we use Fock representation for 
operators a k and at [S2], p3] which associates operator at with some c-number a k and operator 
a k with operator d/da k . Then, with the use of (17.4|) condition (17.7[) takes on the form: 

(E V>VJ- - V Ml w a k , ) A|0) = 0. (7.8) 

\ k r k' / 

It is easy to verify by direct substitution in (17. 8|) that solution of equation (17.81) is written as: 

A 0 = Cexp \ ^ ^2 a t A kk'at > , (7.9) 

W k,k' J 

were C is a constant. To this end it is sufficient to return in (17.91) to quantities a k instead of 
at- Matrix A satisfies the conditions: 

A = M 2 (M*) _1 , A = A T . (7.10) 
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Hence, the ground state energy corresponding to functional A 0 is equal to: 


<0|A+« 0 A„|0) = A£. 


(7.11) 


In Appendix 1 we show that (0 |Aq ifjAolO) = 0. 

From (17.91) . (12.21) . (12.5|) it follows that the wave function of the polaron ground state |^) p 
has the form: 


-jrJ2hka+a k r E fk(a+-a k ) 

= Ce k e k A 0 |0) 


(7.12) 


Accordingly the bipolaron wave function |with regard to (15.3|) . (15.4|) is: 


| i>)bp = Cii}{r)e 


hka+a k R E/fe( a fc 

k (3 k 



(7.13) 


Formula (17. 12ft . (17. 13ft imply that the wave functions of a polaron and bipolaron are delocalized 
over the whole space and cannot be presented as an ansatz (11.21) . 

From formulae (17. 1 2p . (17. 13p it follows that the reason why the attempt of Lee, Low and 
Pines [35j to investigate the polaron ground state energy over the whole range of a variation 
failed was an improper choice of probe function ( 12 .4p which lacks the multiplier corresponding 
to the functional Ao- 

However, it should be stressed that notwithstanding a radical improvement of the wave 
function achieved by introducing the multiplier A 0 in Lee, Low, Pines function enables one to 
take account of both weak and strong coupling, the results obtained by its application are not 
exact. The fact that Tulub function is an ansatz follows from its properties: 


(OlA+ifiAolO) = 0 , E = (0|A+H 0 A 0 |0>, H 0 A 0 |0) = £A 0 |0). (7.14) 


Being an ansatz, Tulubs solution presents a solution of the polaron problem in a specific class 
of functions having the structure of Ao|0). That Tulub ansatz is not an exact solution of the 
problem follows at least from the fact that the use of expression (I3.20|) alone for calculation of 
the energy, for example, in the case of weak coupling yields for E: E = —a — | a 2 |3T] . 

To get an exact coefficient at a 2 in the expansion of the energy in powers of a (11.11) we should 
take into account the contribution of Hamiltonian H \, as the perturbation theory suggests |38j. 

The fact that wave functions (I7.12p . (I7.13|) are delocalized has a lot of important conse¬ 
quences which will be discussed in Section 9. 


20 
























8 Discussion of the completeness of Tulub theory. 


In p3|, |65] a question was raised as to whether Tulub theory [31], [32] is complete. Arguments 
of [S3], |65j are based on the work by Porsch and Roseler eh which reproduces the results of 
Tulub theory. However, in the last section of their paper Porsch and Roseler investigate what 
will happen if the infinite integration limit in Tulub theory changes for a finite limit and then 
passes on to the infinite one. Surprisingly, it was found that in this case in parallel with cutting 
of integration to phonon wave vectors in the functional of the polaron total energy one should 
augment the latter with the addition 5E PR which input will not disappear if the upper limit 
tends to infinity [37], [65]. Relying on this result the authors of [63], [SI] concluded that Tulub 
did not take this addition into account and therefore his theory is incomplete. 

To resolve this paradox let us consider the function A(s) determined by formula (13.1 4 j) 
(accordingly, (13. 19ft in continuous case). As formulae (13.14j) . (13. 19(1 imply, zeros in this function 
contribute into ’’polaron recoil” energy A E given by (j3. 16[) and, according to (j3. 15j) are found 
from the solution of the equation: 




k fk^k 


3 V S ~ U k 


( 8 . 1 ) 


If the cutoff is absent in the sum on the right-hand side of equation (18.Ill , then the solution 
of this equation yields a spectrum of s values determined by frequencies lying between 
neighboring values of and O0k i+1 f° r all the wave vectors /q. These frequencies determine the 
value of the polaron recoil energy: 


AE = \ - u ki) ■ ( 8 - 2 ) 

ki 

Let us see what happens with the contribution of frequencies into A E in the region of the 
wave vectors k where /*. vanishes but nowhere becomes exactly zero. From (18. ip it follows that 
as fk —$■ 0, solutions of equation (18. ip will tend to cu^: zyg —> Accordingly, the contribution 

of the wave vectors region into A E, where fk —y 0, will also tend to zero. 

In particular, if we introduce a certain k° such that in the region k > k° the values of fk 
are small, we will express A E in the form: 


A E 


\ - > 

ki<k° 


(8.3) 
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which does not contain any additional terms. To draw a parallel with Tnlnb approach, there 
we could put the upper limit k° , but no additional terms would appear. 


For example, if in an attempt to investigate the minimum of Tulub functional (13.201) . (15.71) 
we choose the probe function f k not containing a cutoff in the form [ROJJ : 

ft = -Vk exp(-k 2 /2a 2 (k)), 

(8.4) 


a(k ) = 


1 + th 


k h — k 


where a is a parameter of Tulub probe function (15.91) . kt, satisfies the condition a <C kb <C k oc , 
k oc = a\J 3A/4 is the value of the wave vector for which Tulub integral (I5.12|) has a maximum 
[22], |SS3, then with the use of (18.4j) in the limit a —)■ oo , Tulub integral q( 1/A) will be written 
as: 


Q 


5.75 + 6 


exp 


k 2 

K b 


(8.5) 


The second term in the right-hand side of (18.5j) vanishes as k^/a —)■ oo and we get, as we might 
expect, Tulubs result: q( 1/A) ~ 5.75. 

Equation (18. 1 j) . however, has a peculiarity. Even in the case of a continuous spectrum, for 
fk = 0, if k > k° it has an isolated solution u k o which differs from the maximum frequency u: k o 
by a hnite value. This isolated solution leads to an additional contribution into A E: 


- ^fci)+ 5E 


PR 


k{<k° 

SE pr = - W*o) , 


(8.6) 


where z/^o has the meaning of ’’plasma frequency” m- Hence, here a continuous transition 
from the case of f k —> 0 for k > k° to the case of f k = 0 for k > k° is absent. As is shown 
by direct calculation [6?], of the contribution of the term with ’’plasma frequency” 5E PR into 
(18.61) . even for k° —> oo, Porsch and Roseler theory does not transform itself into Tulub theory. 

In Tulub theory we choose such f k which lead to the minimum of the functional of the 
polaron total energy. In particular, the choice of the probe function in the form (18.41) provides 
the absence of a contribution from ’’plasma frequency” into the total energy and in actual 
calculations one can choose a cutoff f k without introducing any additional terms in Tulub 
functional [66|, [67] . 

To sum up, critical remarks in [Blj, [65] are inadequate. Their inadequacy was demonstrated 
in papers }66j . Ea and in work [49] reproduced here. At the present time Tulub theory and 
the results obtained on its basis m , m-m give no rise to doubt. 
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9 Consequences of the existence of translation-invariant 
polarons and bipolarons. 

According to the results obtained, the ground state of a Tl-polaron is a delocalized state of the 
electron-phonon system: the probabilities of electrons occurrence at any point of the space are 
similar. The explicit form of the wave function of the ground state is presented in §7. Both 
the electron density and the amplitudes of phonon modes (corresponding to renormalized by 
interaction frequencies u q .) are delocalized. 

It should be noted that according to (13.151) renormalized phonon frequencies u qi in the case 
of a Tl-polaron have higher energies than non-renormalized frequencies of optical phonons and, 
therefore, higher energies than non-renormalized frequencies of a polaron with spontaneously 
broken symmetry [68] - This holds out a hope to find such phonon modes in experiments on 
light combination scattering and optical absorption. According to [68] . if a polaron (bipolaron) 
is bound on the Coulomb center, i.e. it forms an F-center (F'-center), then all renormalized 
local phonon frequencies u) n have lower energies than the frequency of optical phonon loq does. 
This fact also makes easier experimental validation of the occurrence of delocalized Tl-phonon 
modes with v q . > uj 0 . 

The concept of a polaron potential well (formed by local phonons [68]) in which an electron is 
localized, i.e. the self-trapped state is lacking in the translation-invariant theory. Accordingly, 
the induced polarization charge of the Tl-polaron is equal to zero. Polarons lacking a localized 
’’phonon environment” suggests that its effective mass is not very much different from that of 
an electron. The ground state energy of a Tl-polaron is lower than that of Pekar polaron and 
is given by formula (I5.14p (for Pekar polaron the energy is determined by (j 1.3 [) 7. 

Hence, for zero total momentum of a polaron, there is an energy gap between the Tl-polaron 
state and the Pekar one (i.e. the state with broken translation invariance). The Tl-polaron is a 
structureless particle (the results of investigations of the Pekar polaron structure are summed 
up in [68]). 

According to the translation-invariant polaron theory, the terms ’’large-radius polaron” 
(LRP) and ’’small-radius polaron” (SRP) are relative, since in both cases the electron state 
is delocalized over the crystal. The difference between the LRP and SRP in the translation- 
invariant theory lies in the fact that for the LRP the inequality k c hard < vr is fulfilled, while for 
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the SRP k char a > 7r holds, where a is the lattice constant and k c h ar is a characteristic value of 
the phonon wave vectors making the main contribution into the polaron energy. This statement 
is valid not only for Pekar-Frochlich polaron, but for the whole class of polarons whose coupling 
constant is independent of the electron wave vector, such as Holstein polaron, for example. For 
polarons whose coupling constant depends on the electron wave vector, these criteria may not 
hold (as is the case with Su-Schreiffer-Heeger model, for example |69j) - 

These properties of Tl-polarons determine their physical characteristics which are qualita¬ 
tively different from those of Pekar polarons. When a crystal has minor local disruptions, the 
Tl-polaron remains delocalized. For example, in an ionic crystal containing vacancies, delocal¬ 
ized polaron states will form F-centers only at a certain critical value of the static dielectric 
constant e 0 c- For e 0 > £q c , a crystal will have delocalized Tl-polarons and free vacancies. For 
£ 0 = £q c , a transition from the delocalized state to that localized on vacancies (collapse of the 
wave function) will take place. Such behavior of translation-invariant polarons is qualitatively 
different from that of Pekar polarons which are localized on the vacancies at any value of £$. 
This fact accounts for, in particular, why free Pekar polaron does not demonstrate absorption 
(i.e. structure), since in this case the translation-invariant polaron is realized. Absorption is 
observed only when a bound Pekar polaron, i.e. F-center is formed. These statements are 
also supported by a set of recent papers where Holstein polaron is considered [70]-[72]. The 
approach presented is generalized by the author to the case of Holstein polaron in [75]. 

Notice that the physics of only free strong-coupling polarons needs to be changed. The 
overwhelming majority of results on the physics of strong-coupling polarons has been obtained 
for bound (on vacancies or lattice disruptions) polaron states of Pekar type and do not require 
any revision. 

Taking account of translation invariance in the case of a polaron leads to a minor change 
in the assessment of the ground state, however leads to qualitatively different visions of the 
properties of this state. In paper [32], in the section devoted to scattering of a Tl-polaron, Tulub 
shows that as the constant of electron-phonon coupling increases up to a certain critical value, 
scattering of an electron on optical phonons turns to zero. Hence, when the coupling constants 
exceed a critical value a polaron becomes superconducting. Though in ionic crystals the main 
mechanism of electron scattering is scattering on optical phonons m, it might appear that the 
contribution of acoustic phonons should also be taken into account in this case. However, as 
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it follows from the law of conservation of energy and momentum, a Tl-polaron will scatter on 
acoustic phonons only if its velocity exceeds that of sound 

As distinct from polarons, Tl-bipolarons have much greater binding energy. This leads to 
some important physical consequences. In particular, when a crystal has minor local disrup¬ 
tions, a Tl-bipolaron will be delocalized. Thus, in an ionic crystal with lattice vacancies, forma¬ 
tion of F'-centers by delocalized bipolarons will take place only at a certain critical value of the 
static dielectric constant £q c1 . For e 0 > £ 0ci , the crystal will contain delocalized Tl-bipolarons 
and free vacancies. In the case of £ 0 = £q c1 Tl-bipolarons will pass on from the delocalized state 
to that localized on vacancies, i.e. to F'-center. Such behavior of Tl-bipolarons is qualitatively 
different from the behavior of polarons with spontaneously broken symmetry of Pekar type [8], 
which are localized on the vacancies at any value of £$. 

The fundamental difference between Tl-bipolarons and bipolarons with spontaneously bro¬ 
ken symmetry is that the former are not separable while the latter are separable. This is due to 
the fact that in the case of bipolarons with spontaneously broken symmetry interaction between 
electrons and polarization has the form: <f> (fq, r 2 ) = F(rq) + F(fq)- For |rq — r 2 1 R, where R 
is the bipolaron radius, bipolaron equations separate into two independent polaron equations. 
This fact enables us to treat a bipolaron state with spontaneously broken symmetry as a bound 
state of two polarons p|. In the case of Tl-bipolarons: $(^ 1 ,^ 2 ) = $(rq — r 2 ). I 11 this case, 
splitting of the bipolaron interaction functional into the functionals of interaction of individual 
polarons is impossible for any |fq — r* 2 | and treatment of Tl-bipolarons as composite states be¬ 
comes invalid. This conclusion corresponds to modern ideas that a quantum-mechanical system 
cannot be separated into independent subsystems m 

For zero total momentum of a bipolaron, Tl-bipolarons, being delocalized, will be sepa¬ 
rated from those with broken translation invariance by an energy gap. As with Tl-polarons, 
in the case when the coupling constant exceeds a certain critical value, Tl-bipolarons become 
superconducting. As is known, interpretation of the high-temperature superconductivity rely¬ 
ing on the bipolaron mechanism of Bose-condensation runs into a problem associated with a 
great mass of bipolarons and, consequently low temperature of Bose-condensation. The pos¬ 
sibility of smallness of Tl-bipolarons mass resolves this problem. It should be stressed that 
the above-mentioned properties of translation-invariant bipolarons impart them superconduct¬ 
ing properties even in the absence of Bose-condensation, while the great binding energy of 
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bipolarons substantiates the superconductivity scenario even in badly defect crystals. 

10 Conclusive remarks. 

At the present time Tulub theory and quantitative results obtained on its basis give no rise to 
doubt. The quantum held theory under consideration is nonperturbative and can reproduce 
not only the limits of weak and strong coupling but also the case of intermediate coupling. 

One of the most effective methods for calculating polarons and bipolarons in the case of 
intermediate coupling is integration over trajectories m- Unless properly modified, this ap¬ 
proach is not translation-invariant since in this method the main contribution into the energy 
levels is given by classical solutions (i.e. extrema points of the exponent of classical action, 
involved in the path integral), ffowever, such solutions, in view of translation invariance, are 
not isolated stationary points, but belong to a continuous family of classical solutions obtained 
as a result of action of the translation operator on the initial classical solution. Accordingly, 
the stationary phase approximation is inapplicable in the translation-invariant system. In the 
quantum held theory some approaches are developed for restoring translation invariance. They 
are based on introducing collective coordinates into the functional integral [TFj. However, they 
have not been used in the polaron theory as yet. Therefore it is not surprising that the method 
of integrals over trajectories employed in the plaron theory yields a result coinciding with the 
semiclassical theory of the strong coupling polaron US]. 

Recently in the polaron theory a powerful computational method, namely Monte-Carlo 
technique has been developed [80], [8TJ. This procedure, being only a calculation tool, cannot 
reproduce the results of Tulub ansatz unless properly modified. As for Monte-Carlo diagram 
technique, the obstacle to checking Tulub ansatz in the strong coupling limit is presented by 
the necessity to calculate diagrams of very high order. 

To sum up, Pekar ansatz (II.2p provides an original assumption of the form of the solution 
which was confirmed in the course of numerous examinations. For nearly eighty-year history 
of the polaron theory development (if dating from Landau short paper [82]) ansatz (1 1.2 [) has 
been recognized to be an asymptotically exact solution of the polaron problem in the strong 
coupling limit. 

Tulub ansatz (§7) provides another assumption of the form of the solution whose structure 
is determined by the form of the function Aq| 0). In terms of this assumption Tulub solution is 
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also asymptotically exact. Since Tulub solution yields a lower energy value for a polaron, from 
the variational standpoint, preference should be given to Tulub ansatz. 

Hence, polaron theory in no way can be considered to be complete. In the framework of 
Tulub ansatz great work is to be done to revise many concepts (such as superconductivity) in 
condensed matter physics. Extension of the application area of Tulub ansatz to other divisions 
of the quantum field theory can lead to radical revision of many results which nowadays seem 
doubtless and vice-versa. Thus, for example, non-separability of a bipolaron state in the polaron 
model of quarks [83] (the role of phonons in [S3J is played by a gluon field) provides a natural 
explanation of their confinement. In paper it is noted, that in Tl-theory there is no need 
to use Higgs mechanism of spontaneous symmetry breaking to get the elementary particles 
masses. 

The author expresses sincere gratitude to Prof. A.V. Tulub for numerous discussions and 
valuable advices. The author is also grateful to N. I. Kashirina for valuable discussions. 
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Appendix 1. 


Hamiltonian H\ involved in (13.2p has the form: 
^ ^ (hfc~t~ f fc/ied/c)(h/c 


_^ A: 

k ) + ~~~ fk'( a k a k&k' a t' a k) ' 

k.k' 


1 

2 m 


Wa^a^akak' , (Al.l) 

h h h 1 ~ " h h f 

Ay Ay , Ay Ay , Ay 

where hujk is given by expression (13.51) . Let us apply the operator Hi to functional A 0 (17.91) . 
We will show that (OIA^LRAqIO) = 0. Indeed, the action of A 0 on Hi terms containing an odd 
number of operators in H\ (i.e. the first and second terms in Hi) will always contain an odd 
number of terms and mathematical expectation for these terms will tend to zero. 

Let us consider mathematical expectation for the last term in H \: 


<0|A 0 + V kk'a^a^,a k a k ' A 0 |0). (^l- 2 ) 

k,k> 

The function (0 |Aq a^a^,akak'A. o \0) represents the norm of vector afcafc/A 0 |0) and will be posi¬ 
tively defined for all k and k!. If we replace k —y —k in (Al.2) than the whole expression will 
change the sign and, therefore, (Al.2) is also equal to zero. Hence (0|Ag LRAolO) = 0. 
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Appendix 2. 


Let us show that (15.41) . (15.5)) follow from (15.1)) . (15.2)) . To this end notice that analytical proper¬ 
ties of D(s) pointed out in [ 32 ] straightforwardly follow from (I3.19p . Indeed, the pole of D(s) 
can lie only on the real axis since in view of positive dehnedness of u k k A /| in (13.191) equation: 


i + JL = 

37T- J ( uf. - So) 2 + e 2 


(-42.1) 


can be fulfilled only for e — 0. Besides, D(s ) is a monotonously increasing function s since for 
s < 1, D'(s) > 0, and as So —>■ oo, D(s) turns to unit. Therefore D(s) cannot have zeros for 
—oo < Sq < 1. Hence function (s — l)D(s) can be presented in the form: 


ds' 


(s — l).D(s) 2m J (s'— s)(s'— l).D(s') ’ 

C-\-p 


(42.2) 


where the contour of integration in Caushy integral (.42.2) is shown in Fig.2. The integrand 
function in (42.2) has a pole for s' = 1 and a section from s' = 1 to s' —» oo. Performing 
integration in (42.2) with respect to the upper and bottom sides of the crosscut we will get 
integral equation (15.51) . 

Appendix 3. 


Let us perform a detailed calculation of the quantity A E (I5.7[) with the use of probe function 

m>- 

To this end, to calculate the real part of D{ujp) involved in (15.71) we use Sokhotskys formula: 


= v 7jr—2 + > 

Uk ^p ^ ^k k'p 


ReD(c<2 2 ) — 1 + 


P 2,4^ u k 




v , - ■ 37r 2 J ' u 2 _ u 2 


dk . 


It is convenient to present Re D in the form: 


Re D — 1 T I\ T I 2 , 


h = 


2iA 


37T 2 


m 


dk 


Uk + U p 


, h = v 


(jJ<r 


flk A dk 


37T 2 J (u k - uj p ) ( u k + Up) 

0 
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Substituting f k in the form of (15.91) into these expressions we present Ji as: 


h — 


8a 


3V2 


7T 


■ e -^ dk .^±3 f dk 

3 V^rt J k 2 + p 2 + 4 


Assuming k/a = k in the strong coupling limit (a —> 00 ) we write for ip 


/1 = 


8aa 


3V2vr 

Accordingly, J 2 takes on the form: 


71" 7T _ -2 


— - -pe p | 1 - — / e df 


7T 


/, = p 


e k2 / a2 k 2 dk 


This integral can be expressed as: 


37T\/2 i ^p){p^k + 


A — Ao + Ai j 


where: 


lQau p 
An 


I 


37T\/2 V P 2 + 2 / J k 2 + p 2 + 4 


0 —k 2 /a? 


-dk , 


Ai — 


00 , o , . 

16an; p (a; p -l) p /' e _fc /“' 


3nV2(p 2 + 2) J k 2 -p 2 
The integrals involved in A 0 and J 2 i will be: 


dk. 


-k 2 /a 2 1 

-dfc = 


k 2 + p 2 + 4 a 


e dt 


7T 


7T 


2 p ’ 


A 


g— k 2 /a 2 

k 2 — p r ‘ 


-dk = 


n e 


-p 


a p 


e l dt 


As a result, J 2 has the form: 


A = 


2 aap ^2 


3 y/2 

Finally, ReP will be written as: 

Re D = 1 + 


e dt 


7T 


f e^dt . 
3\/2vr ./ 


4aa 


3\/27r 


1 — pe p e dt — pe p e dt 
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This result reproduces the quantity given by formula 05.101) . For the imaginary part Im D by 
Sokhotskys formula, we get: 


Im D = 


f%k*ui k 6(u>% - ujp)dk = -^flp 3 . 


3n 


P ' 67T JP1 


As a result, \D(oj 2 )\ is expressed as: 


\D\ 2 = (Re D) 2 + (Im D) 2 = £ a 2 a 2 

9 


e P P + — | 1 - p I e~* dt- pe~ p / e* dt 


The first term in formula 05.71) is easily calculated to be: 


i r k fk „ 3 2 

4n 2 J (1 + Q) 16 


In calculating the second term in 05.71) we will separate out integral I p in it: 

-k 2 /a 2 + OJkipJk + ^p) + 1) 


I - / e -fc a /° — 


(uk + UJp) 2 


As a —y oo, it is equal to: 


! p -a 3 4 


7T 


1 - p 3 e p / e ^dt (2 + 4p 2 ) + 2p 4 1 = j " (1 - 0(p)), 


a 3 x /7T. 


where: 


fi(p) = 2p < (1 + 2 p 2 )pe p2 / e * 2 dt — p 2 


which corresponds to the expression h2(p) in 05.121) As a result, the second term in formula 05.71) 
will be: 


1 47ra 
12tt 4 ^2 


4 


ipP7 3 


UJn 




dp. 


As a —y oo, this expression takes on the form: 


a 2 a 4 


371^ , 


(i - n(p)) 


3 -p 2 


I^K)I 2 


3 

16 


dp = —a 2 <? , 


where q = g(0) is given by expression 05.12|) . Hence, finally for A E (I5.7|) we get: 

AE = Jd a2 ( 1 + q ') ’ 

16 

which corresponds to the first term in the right-hand side of 05.lip . 
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